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1 Introduction 



Abstract 

Coined quantum walks may be interpreted as the motion in position space of a quantum particle 
with a spin degree of freedom; the dynamics are determined by iterating a unitary transformation 
which is the product of a spin transformation and a translation conditional on the spin state. 
Coined quantum walks on 1 d can be treated as special cases of coined quantum walks on W l . We 
\& \ study quantum walks on M. d and prove that the sequence of rescaled probability distributions in 

position space associated to the unitary evolution of the particle converges to a limit distribution. 

o 

S3 

P^ 

Several kinds of quantum walks have recently been studied by many authors; a nice overview is given 
by JO]. So-called "coined quantum walks" on finite graphs, introduced in pQ, are proving to be of 
some interest in quantum informatics, where they have been used to devise fast search algorithms 
Pl llOlUB] , Quantum walks of the type considered in this article, namely, coined quantum walks on 
Z d or M. d , were first introduced in jj, though |2] and ^Sj can be considered as precursors in some 
respects. 

A coined quantum walk describes the evolution of a quantum system under iteration of a certain 
kind of unitary map. The state of the system is a vector in a product Hilbert space, having position 
degrees of freedom and an internal degree of freedom such as spin or polarization. A step of the 
walk consists of a unitary transformation of the spin degree of freedom, which one may think of as 
"tossing a quantum coin" , followed by a translation conditional upon the state of the coin. Here 
we restrict our attention to walks on Z d or M. d and we prove that the position distribution of the 
quantum walker, properly rescaled, converges as the number of steps tends to infinity. This type 
of convergence was first discovered by N. Konno |lll H^] in the case of one-dimensional lattices Z. 
A different proof of Konno's theorem was given in |7j, generalizing the convergence to quantum 
walks on higher dimensional lattices Z d . In this article we present a further generalization of the 
result to quantum walks on R d , and we obviate a technical difficulty of 7 that required an extra 
hypothesis there. 
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2 Basic definitions and review of Konno's theorem 

We start by reviewing some basic definitions and stating Konno's theorem for a coined quantum 
walk on Z. 

A simple quantum walk on Z can be denned as a sequence 

tl>, u*/>, U 2 il>, c/ 3 v, • ■ ■ , (1) 

of unit vectors in the Hilbert space £ 2 (Z) <g> C 2 obtained by iterating a unitary operator U of the 
form 

U = S(I®C). (2) 

In (J2J), J denotes the identity operator on £ 2 (Z), C denotes a unitary operator on C 2 , and S denotes 
the "conditional shift" operator 

S(s m ® ei) = s m+ i ® ei 

S(s m (g) e 2 ) = s m _i <g> e 2 , (3) 

where e\ and e 2 are the standard basis vectors for C 2 and s m is the vector in £ 2 (Z) whose m th 
member is 1 and all other members are 0. Since the vectors of the sequence Q are normalized, the 
numbers 

P n (m) = |(s TO (g)ei, C/"V)| 2 + | (s m ®e 2 , U n ip}\ 2 (4) 

satisfy J2 m Pnijn) = 1 and thus define a sequence of probability measures on Z. Note that the only 
variable considered here is the position of the walker; the internal degree of freedom is disregarded 
by summing over both states in the coin space. Konno's theorem |12| states that the probability 
measures 

P n (m) 5 m/n (5) 
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converge weakly as n — > oo to a probability measure that depends on the initial state ip (here 
6 m / n denotes a point-mass at m/n). 

The dynamics described in (Q, @, (j2Jl is reminiscent of simple random walk on Z, though 
there are a few important differences. As in ordinary random walks, the rule @ for stepping left 
or right is the same at all locations m 6 Z, that is, the process is spatially homogeneous. Unlike 
ordinary random walks, the probability distributions (@J) for the system's position are not related 
by Markov transitions. Similarly to random walk, the sequence of rescaled position distributions 
l|^)l converges weakly, but unlike random walks, the weak limit is obtained by rescaling by a factor 
of 1/n rather than l/\/n and the limit distribution depends on the initial state (and is not a normal 
distribution). 

In the next sections we state and prove a generalization of Konno's theorem in which (i) the 
walk takes place in <i-dimensional space M. d , (ii) any finite number s of conditional translations 
are allowed at each step, and (iii) the translations are not assumed to generate a lattice in R d . 
Quantum walks on the ci-dimensional lattice Z d may be viewed as special cases of quantum walks 
on R d where the shifts generate a discrete lattice. 

3 General quantum walks convergence theorem 

Broadly speaking, a quantum walk is characterized by the unitary operator U that generates the 
walk, U being a particular kind of unitary operator on a Hilbert space L 2 (Z d ) <g> C s or L 2 (R d ) (g) C s , 



the former for quantum walks on Z d and the latter for quantum walks on R d . We now focus our 
discussion onto quantum walks on M. d , for quantum walks on 1 d are easily embedded into the 
framework of walks on M. d . 

Let U = L 2 (R d ) (g) C s . There is a natural isomorphism between L 2 (R d ) ® C s and L 2 (R d ,C s ); 
it will be helpful to represent members of Ti as members of L 2 (M. d ,C s ), and indeed to represent 
members of C s as column vectors for the purpose of interpreting matrix operations. Let C denote 
a unitary "coin tossing operator" on the "coin space" C s and let I s denote the identity operator on 
C s . Let ei, e2, ■ • • , e s be the standard ordered basis of C s . The steps of the quantum walk involve 
a "conditional translation" operator T on Ti defined by 

T(ip®ej) = T v ^®e 3 , (6) 

where T v . denotes translation by Vj in Z d or R d , i.e., (T v .ip)(x) = ip(x — Vj). The vectors Vj are 
arbitrary fixed vectors in Z d or M. d . A step of the quantum walk is effected by a unitary operator 

U = T(I®C), (7) 

where / denotes the identity operator on L 2 (Z d ) or L 2 (R. d ). 

The operator U generates the quantum walk. In the Schrodinger picture, the quantum walk 
may be thought of as a sequence 

p, UpU*, U 2 pU* 2 , U 3 pU* 3 , ... (8) 

of density operators on Ti. (Recall that a density operator is a nonnegative trace class operator with 
normalization Tr(p) = 1 |g). The density operator p for the normalized pure state ip £ Ti has integral 
kernel ip(x)ip(y). Thus, the definition (jHJ) of quantum walk in terms of density operators generalizes 
the original definition QJ for pure states ip.) The probability measures J2J can be expressed in 
terms of density operators instead of wavefunctions as follows. Given an initial density operator p 
on Ti and any n £ N, there is a probability measure P p ,„ on Z d or M. d that may be defined as a 
linear functional on Co(Z d ) or Co(K d ) by the formula 

P p ,„/ = Tr(U n pU* n (m[f]®I s )), (9) 

where Cg denotes the set of continuous functions that tend to at oo, / is an element of this set, 
and m[f] is the operator of multiplication by /. Formula ^ may be written as 

P p .nf = Tr(Tro (U n pU* n )m[f}), (10) 

where Trc= denotes the partial trace; for the special case studied in Section [21 one may verify 
that P Pi „ is the same as P„ of |@J when p is a pure state ip- The quantum walk as such is the se- 
quence of density operators © and the corresponding probabilities (|10l) , which express probabilitcs 
concerning the position of a "particle" without regard to its internal "spin" state in C s . 

One observes that the probability measures P Pi „ spread out in Z d or M. d as n increases. However, 
we will discover that dilating each of the P Pi „ by a factor of 1/n produces in a sequence of probability 
measures that converges weakly. In other words, if X n is a random vector describing the position 
of the quantum walk after n steps, then the sequence X n /n converges in distribution. 

Now we are prepared to state a generalized version of Konno's theorem ^] El E In the 
following, Cb(R d ) denotes the space of bounded continuous functions on M d , and a sequence of 
probability measures is said to converge weakly if the measures converge pointwise as functionals 
on C b (R d ) 0. 



Theorem 1. Let U be a unitary operator on Ti — L 2 (R d ) ® C s of the form (J7J). 

There exists a spectral measure P{dx) on R d such that for any density operator p on Ti, the 
probability measures ¥ p "n defined by 

P$/(*) - Tr(Tr Cs (U n pU* n )rn{f(x/n)}) (11) 

converge weakly as n — > oo to t/ie probability measure W P (E) = Tr(pP(E)). 

We will prove this theorem — and somewhat more — in the next section, but first it is ap- 
propriate to explain in which sense Konno's theorem for walks on the lattice is generalized by 
Theorem ^ Convergence of quantum walks on M. d implies convergence of quantum walks on Z d 
because quantum walks on Z d can be identified with quantum walks on R d that have integral steps. 
Indeed, the space L 2 (Z d ) can be identified with the subspace K C L 2 (M. d ) consisting of functions 
that are constant on each unit cube centered at a lattice point. The subspace K. is invariant under 
translation by any lattice vector in Z d , so that a quantum walk on Z d may be viewed as a quantum 
walk on K d whose initial density operator is supported on the subspace K(g> C s of L 2 (M. d ) (g> C s . The 
analog of an initial state localized at a vertex of the lattice is a state supported on a cube centered 
at the corresponding lattice point. The limit distribution for the discrete case can be retrieved 
by taking the initial state in M. d to have support on the unit cube centered at the corresponding 
lattice point. Convergence of the P p ™n for the walk on R d then implies the convergence considered 
in [T21I7| for walks on Z d . 

Remark 1. The scaling by 1/n relates the description of the asymptotic distribution of the quan- 
tum walk to the equation of motion for a quantum particle whose wavefunction propagates linearly 
in time, with a constant (but random) velocity vector that has the distribution P p . Hence, a phys- 
ical interpretation of the commuting bounded operators V\ , . . . , Vj, corresponding to the spectral 
measure P (cf. Theorem EJl is that they are commuting observables giving the asymptotic velocity 
vector of the quantum random walk. Note, however, that although this describes the asymptotic 
distribution, it does not describe the quantum random walk precisely; indeed, the observables 
corresponding to the positions at different times do not commute and hence one cannot regard 
the quantum random walk as a stochastic process, since that would mean one could introduce a 
common probability space for the positions at different times. 

4 Proof of Theorem [T] 

Theorem n is formulated in the Schrodingcr picture, where the evolution is applied to the state of 
the system. In the Heisenberg picture, where the the state of the system remains constant and the 
physical observables evolve, the quantum walk concerns the map 

X i — > X, U*XU, U* 2 XU 2 , U* 3 XU 3 , ... 

from bounded operators X on Ti to sequences of bounded operators. We are interested in operators 
X that represent physical observables of position alone. These are operators of the form m[f(x)]®I s , 
where m[f(x)] denotes the operator ip(x) i — ► f(x)ip(x) on L 2 (M. d ). Theorem-His the consequence 
of a stronger, dual formulation in the Heisenberg picture. In the following, the notation s • lim 
designates the limit in the strong operator topology, the topology in which nets of operators converge 
if and only if they converge pointwisc as functions on Ti. 

Theorem 2. Let U be a unitary operator on Ti — L 2 (R d ) ® C s of the form (J7J). 



There exists a spectral measure P[dx) on M. d and a corresponding family of commuting bounded 
self-adjoint operators V\, . . . , Vd such that 



s-]imU* n (m[f(x/n)]®I s )U n = / f(x)P(dx) = f{V u ...,V d ). 

n->oo J Rd 

for all f e C b {R d ). 

Theorem n follows from Theorem [5] and the observation that if p is a pure state, then 

Urn P$J(x) = lim Tr(U n pU m {m[f(x/n)]®I s )) 

= Tr (p (s • lim U* n (m[f(x/n)} <g> J S )C/™) 



(12) 



Tr(p ff(x)P(dx)) = f f(x)Tr( P P(dx)) 



(13) 



for any / E Cb(M. n ). The validity extends to general p by a density argument. 

Thus, it suffices to prove Theorem [21 Here is the plan of the proof: Note that, for each n, the 
map 

U n :w^ U* n (m[e lx - uj/n ] <g> I s )t/ n (14) 

is a strongly continuous unitary representation of K d on L 2 (R d ,C s ), i.e., U n is a group homomor- 
phism from R d to the unitary operators on H such that XJ n (tu)ip is continuous in w for each tp e W. 
We will prove that the U n (w) converge pointwise in the strong operator topology as n — > oo 
and we will identify the limit as a uniformly continuous unitary representation of M d . By Stone's 
spectral theorem there exists a spectral measure P(dx) on M. d such that 



s • lim U* n (m[e iu - x/n ] <8> I s )U n = 

n — >oc 

for all w € R d . This proves that 

s-lim[/*™(™[/(x/n)](g)/ s )t/ 11 = 



e luJ - x P(dx) 



f(x)P(dx) 



(15) 



for all / G span{e lw:l:: : to e K d }. Finally, to complete the proof of Theorem|21 we will extend the 
convergence in (|15f) to all functions / G Cb(M. d ). 

We begin by showing that U„(oj) converges strongly at each u> £ M. d . To do this, we will first 
prove that the infinitesimal generators of these unitary groups converge on a dense subset of Ti to 
a bounded skew-Hermitian operator, and then invoke the Trotter-Kato Theorem. 

For arbitrary but fixed lo £ M. d , the one-parameter unitary group {XJ n (tui)}t^s, has infinitesimal 
generator 

G n = \U* n {m[iu ■ x] (8 I s )U n . (16) 

Observe that if D{u) is the operator on C s represented by the diagonal matrix 



D(w) 



uj ■ v\ 

UJ ■ V2 



LO ■ V s 



then, for any / G L 2 (R d ), 

T*(m[uj ■ x] <g> I 8 )T(f (gi e,-) = m[w ■ (a; + «,)] ® 4(/ ® e,-) = f m[w ■ x] <8> I s + 1 <g> I>(w)J (/ <8> e,-) 
and hence 

[/'(m[iwi]8/ s )t/ = (I<g)C*)T*(m[iw-a;]® J s )T(J<g)C) 

= m[iw-x]®I s + i(I®C*D(ui)C) . (17) 

Applying (|17|l recursively in (|16fl shows that 

_. n — 1 

£ n = imfiwalg/, + * -y^U* j (I®C*D(uj)C)U j . (18) 

n z — ' 

3=0 

The first term on the right-hand side of (|18fl is an unbounded multiplication operator; as n — > oo 
these operators converge to the zero operator on the dense subset of L 2 (R d ,C s ) consisting of 
functions of bounded support. The second term on the right-hand side of (|18(l is a bounded 
operator; we will prove that these bounded operators converge strongly as n — ► oo. 
We take Fourier transforms to identify the limit of the operators 



n-l 



i ■- j - 
-J2u* j (I®C*D{u))C)U j 



(19) 



Let T denote the Fourier transform 



on W l . This is a unitary transformation from L 2 (M. d ) onto an isomorphic space L 2 (R d ) with inverse 

T*(<j>)(x) = (2ixY d l 2 f 4>{k)e lk - x dk . 



The spaces L 2 (R d ) ® C s and L 2 (R d ) ® C s are mapped onto one another via the unitary transfor- 
mations T <£> I s and T* I s . We will denote T ' ® I s by T from now on. The operator U of is 
unitarily equivalent to the operator U = TUT* on Ti. We will think of Ti (the momentum space 
of the system) as L 2 (M. d , C s ), that is, we represent vectors in Ti by square-integrable column- vector 
valued functions on M. d . Then the operator U may then be represented by a "matrix- multiplication 
operator" 

(U<f>)(k) = U k <P(k) 

with 



U k = 



-ik-vi 



o 





-ik-V2 






-ik-Vs 



C\\ C12 

C21 C22 

Csl C S 2 



C*2s 



(20) 



■■■ e 

where the Cjk denote the matrix elements of the coin toss operator C . The Fourier transform of 
the operator — 5Z"T J7* J '(7® (C*D(a;)C))t/ : ' is the matrix-multiplication operator 



ra-1 



0(fc) 



1 ^ ^ 



(21) 



i=o 



For fixed k £ M. d the map X i — ► U^XUk is a unitary operator on the space of s x s matrices 
endowed with the inner product (Y, X) = Tr(Y*X). The Mean Ergodic Theorem implies that the 
matrices in (|21[1 converge as n — ► oo to the orthogonal projection of C*D(u>)C onto the subspacc 
of s x s matrices that commute with Uk- That is, 

-. n—X 

W„{u) = 1™-E^( C *%) C )^ ( 22 ) 

is the unique matrix that commutes with Uk and satisfies 

Ti(N*W k (uj)) = Tr(N*C*D(u>)C) (23) 

for all matrices N that commute with Uk- Since C*D(lo)C is Hermitian and since N commutes with 
Uk if and only if N* commutes with Uk, taking adjoints in (|23[) shows that the matrices Wk(ui) are 
Hermitian. From (|20|l and 122[) it is seen that Wk(co) is a limit of continuous functions of k, hence 
it is a measurable function of k. Moreover, the operator norm of each matrix on the right-hand side 
of (1221) does not exceed the operator norm of D(ui), which is independent of k, and the limit Wk(u>) 
enjoys the same uniform bound. Thanks to these bounds, Lebesgue's Dominated Convergence 
Theorem implies that the operators (|21|) converge strongly to the matrix-multiplication operator 
W(u>) defined by 

(W(w)0(fc) - W k (w)<t>(k) (24) 

on 7i, a bounded Hermitian operator. Transforming back to operators on Ti, we conclude that the 
operators fH?j) converge strongly to 

V(w) = T*W{lj)T. (25) 

We have now shown that the generators (fT5|) converge to iV(u>) on a dense subset of H. By 
the Trotter-Kato Theorem ^] Theorem 4.5], the operators \5 n {tuj) converge strongly to e ,iV ^- ) 
for each tgl. In particular, U n (w) converges strongly to e lVlyUJ \ i.e., 

s-limU* n {m[e ZUJ - x/n ] ® I s )U n = e lV{uj) . 

n — >oo 

It follows from l|22() . I|24() and (|25|l that V{u>) depends linearly on w e l d ; in particular, uj \ — > 
V(u>) is continuous. Moreover, lo i — ► e lV (") is a d-parameter group of unitary operators on Ti 
because each u> i — ► U„(w) is. It follows that e'^^' is a uniformly continuous d-parameter unitary 
group. Stone's spectral theorem for one-parameter unitary groups can be generalized to d-parameter 
unitary groups: for any strongly continuous d-paramcter unitary group U(w) there exists a spectral 
measure P(dx) on M. d such that 

U(w) = / e luJ - x P(dx) 

JR d 

(see, e.g., Chapter 2.4]). For one-parameter groups e iu ' v , P(dx) is just the spectral measure 
associated to the self-adjoint operator V, and in general P is the spectral measure of a d-tuple 
(Vi, . . . , Vd) of commuting bounded operators. Letting P(dx) denote the spectral measure associ- 
ated to e lV ^\ what we have shown up to this point is that 

s-\imU* n {m[e luJ - x/n ]®I s )U n = e iV(w) = / e iuJ - x P(dx) (26) 

n->oo J Rd 

for all u> £ R d . 



To complete the proof of Theorem[21 we have to extend to all functions in Cb(R d ) the convergence 
H2tj|) just established for functions e luJ ' x . By linearity, the assertion of the theorem holds for all 
/ G span{e Ma: : co G R d }. Let / be an arbitrary but fixed function in Cfc(M d ) and define 

L[/,n](V) = U* n (m[f(x/n)]®I s )U n ij 

on Ti. We want to prove that 

lim L[f,n](ip) = [ f(x)P(dx)ip 

for all ip ^7i. The sequence {L[/, n]} is equicontinuous on Ti thanks to the uniform bound 

\\U* n {m[f{x/n)]®I s )U n \\ = H/IU, 

so it suffices to prove the convergence of the sequence {L[/, n]} on a dense subset of Ti. We will 
show that h[f,n](tp) converges for all ifi that have bounded support. 

We continue to represent members of Ti by functions in L 2 (R d ,C s ). Let tp G Ti have bounded 
support, in the sense that there exists r > such that ip(x) equals G C s when |x| > r. We claim 
that h[f,n]ip tends to a limit. We may assume that ip has norm 1 without loss of generality. Let 

v = max {|f 7 |} , 

l<j<s u ' 

where the Vj are the translation vectors appearing in the definition |jS| of T. Then U n tp is supported 
on the ball of radius r + nv in R d . For any g G Cb(M. d ), 

\\L[f,n]i>-L[g,n]iP\\ = \\U* n ((m[f(x/n)} - m[g{x/n)]) ® I s )U n ^\\ 

< sup {|/(aO-ff(a:)|} (27) 

\x\<r-\-v 

for all n because of the way that / and g are scaled. By the Stone- Weierstrass Theorem, the linear 
span of {e luJ ' x : u> G M. d } is dense in the space of continuous functions on any compact subset of 
M. d , so f{x) can be uniformly approximated within arbitrary e > on the ball of radius r + v 
by some function g e G spanje 4 "' 1 : to G R d }. We know that each sequence L[g e ,n]ijj converges as 
n — > oo. On the other hand, L\c/ e ,n]ip is within e of ~L[f,n]ip uniformly in n by l|27|) . It follows 
that {L[f,n]'ip} converges as n — > co. The limit is a continuous function of / and hence must 
agree everywhere with J f(x)P(dx)ip = f(Vi, . . . , Vj). This proves (|12() for arbitrary / G Ch(R d ) 
and concludes the proof of Theorem [21 

5 The limit distribution 

In special cases there is a nice characterization of the limit distribution P p of Theorem ^ m terms 
of the distribution of a certain random element of O = M. d x {1, 2, . . . , s}. 

Let p = T pT* be the density operator on Ti corresponding to p. This can be regarded as an 
integral operator p f = L d p (k, k')f(k') dk' with an s x s matrix valued kernel p (k, k'). Since p is 
trace class, the diagonal p (k,k) is well-defined a.e. and integrable. By Theorem^and (|13fl . (|26|l . 
we have l|25|) . 

T? p [e iu - X ] = Tr(p f e %UJX dP{x)\ = Tr(/7e iV(w) ) = Tr(pe tW< - u A = f Tr [e tWk( -^p (k, k)] dk. 



This shows that the characteristic function of the limit P p in Theorem ^ is 

Tr[e tWkiuj) p (k,k)]dk, (28) 



with W k given by (J22J) and lj2lfl) . 

We can better identify the limit measure ¥ p in cases where the matrices U k have s distinct 
eigenvalues Xi(k), . . . , X s (k) for a.e. k S M d . At a point k where there are s distinct eigenvalues, 
let ipj{k) be a normalized eigenvector of U k with eigenvalue Aj(fc), and let P k j = i)j(k)ipj(k)* be 
the corresponding orthogonal projection onto the eigenspace. We may assume that the eigenvalues 
are numbered such that Xj are continuous and differentiable in a neighborhood of k for each j 
0PI-5.4], and we then define 

n(k,j) = i\j(k)V\j(k), (29) 



where V is the gradient with respect to k £ R d . Since |Aj(fc)| = 1, the functions iXj(k)WXj(k) are 
real-valued, and it is an a.e. defined map from M. d x {1, 2, . . . , s} to K d . 

Lemma 1. Suppose that U k has s distinct eigenvalues at almost every k. Then. 

s 

W k (u>) = ^2(u>-n(k,j))P kj (30) 

and thus 

P P [e^] = / y2 eiU ' T(k,J)Tr lP (k,k)Pkj]dk . (31) 

J = l 

Proof. We may assume that the eigenvectors ipj are differentiable in a neighborhood of k |5][II-5.4]. 
Taking the directional derivative d u — us ■ V of U k tpj(k) = Xj(k)ipj(k) we obtain, since (|20|l yields 
d u U k = -iD(u})U k , 

-iD(u)Uk1>j(k) + U k d^(k) = (d^X 3 (k))^(k) + X 3 {k)d^ 3 {k). (32) 

Further, 

^(k)*U k = (UZ1>MY = (XjW^ik))* = X 3 (k)^(kr 
and thus P2"|l implies 

^•(fe)*(-i£»(o;)C/ fe )^-(fc) = Vi(fc)*(^A,-(fc))^(fc) = 9<A-(A;). (33) 

Since Wfc(w) commutes with [/& it is of the form V ■ ajP k j. Taking N = N* = ipj(k)ipj(k)* in J33J) 
yields, using C*D{uj)C = U^D(uj)U k , jS3> and (EHJ, 



Tr 



^•(fc)^i(fc)*Wib(w)) = Ir(^(fe)^(fe)*C*£>(a;)c) = ^{kfU^D^U^^k) 



= Xjik^jik)* D{u)U k il) {k) = iXj(k) d„\j(k) = w • 7r(fc, j). 
This gives (J3D|, and (EU then follows from J2jj&. D 



Let P p denote the probability measure 

F p (dk,j) = Tv[p(k 1 k)P k ,]dk (34) 

on £1 = M. d x {1, 2, . . . , s}. Then the right-hand side of H31|) is the characteristic function of the 
induced probability measure P p o 7r _1 on R d , which identifies this probability measure as P p . In 
other words, if Y is a random element of £7 with the distribution P p , then P p is the distribution of 
the random variable ir(Y). We conclude with an interesting special case: 

Corollary 1. Suppose that the initial state is is of the form p ~ po -I s , a tensor product of a 
position density operator po on L 2 (R d ) with a maximally mixed coin state. Suppose further that the 
matrices Uk have s distinct eigenvalues Ai(fc), . . . , X s (k) almost everywhere on M. d . Then the limit 
P p is the distribution of the random variable tt(Yq,Z), where it is given by (I29|) . Yq is a random 
vector in M. d with density function Trpo(fc, k) and Z is uniformly distributed on {1, . . . , s}, with Yq 
and Z independent. 

Proof. We have p = p~o ® -I s and, for the kernel, p (k, k) = p~o(k, k)-I s . Since each Pkj has rank 
1, 1(3*4^1 shows that 

V p (dk,j) = ±p a (k,k)dk. 

This is a product measure and thus, if Y above is written as (lo, Z), then Yq and Z are independent 
with the stated distributions. □ 

Note that in Corollary ^ the definition of the quantum walk, as given by JSJ and (JJJ) , affects 
only 7r, while the initial state affects only the distribution Yq. 

Remark 2. As explained at the end of Section|3 the results transfer to quantum random walks on 
Z d too. It can be verified, using the method described there, that in this case the limit distribution 
is described by the formulas above, but with k 6 K d , the dual group (a d-dimensional torus). This 
generalises results in [2j- In particular, if the initial position is G II 1 with a maximally mixed coin 
state, Corollary Q] holds with Yq uniformly distributed on K d . 
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